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The conservative dynamics of gravitationally interacting two-point-mass systems has been re¬ 
cently determined at the fourth post-Newtonian (4PN) approximation [T. Damour, P. Jaranowski, 
and G. Schafer, Phys. Rev. D 89, 064058 (2014)], and found to be nonlocal-in-time. We show 
how to transcribe this dynamics within the effective one-body (EOB) formalism. To achieve this 
EOB transcription, we develop a new strategy involving the (infrnite-)order-reduction of a nonlocal 
dynamics to an ordinary action-angle Hamiltonian. Our final, equivalent EOB dynamics comprises 
two (local) radial potentials, A(r) and D{r), and a nongeodesic mass-shell contribution Q{r,pr) 
given by an infinite series of even powers of the radial momentum pr. Using an effective action 
technique, we complete our 4PN-level results by deriving two different, higher-order conservative 
contributions linked to tail-transported hereditary effects: the 5PN-level EOB logarithmic terms, 
as well as the 5.5PN-level, half-integral terms. We compare our improved analytical knowledge to 
previous, numerical gravitational-self-force computation of precession effects. 
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I. INTRODUCTION 

The impending prospect of detecting the gravitational- 
wave signals emitted by coalescing binary systems gives a 
new incentive for improving our theoretical knowledge of 
the dynamics of two-body systems in general relativity. 
Recent developments have shown that a useful strategy 
for accurately describing the dynamics of binary systems 
is to combine, in a synergetic manner, the information 
gathered from several different approximation methods: 
notably, the post-Newtonian (PN) formalism, the grav¬ 
itational self-force formalism, full numerical relativity 
simulations, and the effective one-body (EOB) formal¬ 
ism. In a recent paper [1] , we have succeeded in deriving 
the conservative dynamics of a two-body system at the 
4th post-Newtonian (4PN) approximation [i.e., including 

fractional corrections of order (v^/c^+GM/{c^r)'j to the 
Newtonian dynamics]. (Our work was the culmination of 
previous partial, 4PN-level, investigations m!.) 

The aim of the present work, is to transcribe the 
Taylor-expanded 4PN dynamics of Ref. [T] within the 
EOB formalism [TOUI^ . The EOB formalism provides an 
analytical framework for the description of the relativis¬ 
tic two-body problem which has many useful features: 
notably, (i) it encompasses natural resummation tech¬ 
niques allowing one to extend the validity of perturbation 
results up to merger; (ii) it can extract nonperturbative 
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information contained in a few numerical relativity sim¬ 
ulations with analytic perturbative information; and (iii) 
it provides accurate gravitational waveforms correspond¬ 
ing to the full coalescence process from early inspiral to 
ringdown. (For a sample of recent EOB results see Refs. 
[I4lll8j .) Until now, the EOB description of the two- 
body dynamics is fully known only at the 3PN level [12] , 
though some parts of the EOB description are known to 
higher PN accuracies. For instance, the main EOB radial 
potential A(r) is analytically known, to linear order in 
the mass ratio, up to the 9.5PN level HHID], while the 
logarithmic contributions to the secondary EOB radial 
potential D{r) are known up to the 5PN level 


II. EOB REMINDER 


To set the stage, let us briefly recall the main features 
of the EOB formalism. The basic idea is to map the 
conservative relative dynamics of a binary system (seen in 
the center of mass frame) onto the (equivalent) dynamics 
of an effective body moving in some effective metric g®® 
(with additional, nongeodesic, Finsler-type corrections). 
This is a general relativistic generalization of the well- 
known fact that the Newtonian relative motion of a two- 
body system is equivalent to the motion of a particle of 
mass p := mi m 2 /(mi -I- m 2 ) in the two-body potential 
V{r). Here, mi and m 2 denote the masses of the binary 
system. In the following, we shall also denote 


M := mi + m 2 . 


v := 


M 


mi m2 
(mi +1712)'^ 


( 2 . 1 ) 
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The full Hamiltonian of the two-body system, in the cen¬ 
ter of mass frame (pi -I- P 2 = 0, X 12 := xi — X 2 ), can be 
written as 

7J[xi2,pi;TOi,m2] = Mc^ -|-/r i?"''[xi2, pi; wi, m 2 ], 

J2.2) 

where the reduced, “nonrelativistic” Hamiltonian iJ"’' = 
= {H — Mc^)/fx has a PN expansion of the type 


i7“[Xi2,Pl] = 


1 Pi 


GM 1 


1 


U-nr I ^ trnr 

, I ^^2"iPN + ^4 "2PN 


^ ^3PN + ^ ^4PN + ^ (^ ) • 


We wish to map the “real” Hamiltonian (2.2 1 onto some 


“effective” Hamiltonian i^eff[r^P^] describing a relativis¬ 
tic dynamics, with a general mass-shell condition of the 
type 


+ Qix',P) = o> 


(2.4) 


where 


Qix\p') = 


Q\ 






(2.5) 


denotes contributions which are at least quartic in mo¬ 
menta. As argued in [12] at the 3PN level, and as we 
shall show below at the 4PN level, one can reduce the 
p'-dependence of Q to a dependence on the sole rad ial 
momentum p(.. Then, the mass-shell condition (2.41 is 
quadratic in the time component pg, and iJeff [I'^ P 4 de¬ 
notes the positive-root solution for the effective energy 
^eff — Po • 

At the 4PN level, the a priori unknown functions 
parametrizing the EOB dynamics comprise: 

(i) the two functions A(r'), B{r') parametriz¬ 

ing a generic spherically symmetric metric (in 
Schwarzschild-type coordinates) 

gf^{x')Ax'^^dx''' = -A(r')c2dt'2 

-I- B{r') dr'^ -|- r'^(d0'^ -I- sin^ 9'Aip'^) (2.6) 

[henceforth, we work with the function D{r') := 
{A{r')B{B))-^ instead of H(r')]; 

(ii) the functions ... parametrizing the 

various contributions to Q, Eq. (2.51; and 


(iii) an energy-map function = /(H) relating the 
effective energy = HeB to the real energy 
^reai = H oi the two-body system. 

Following [iniiia, we a priori allow for a general en¬ 
ergy map of the type 


unr / CJnr 

£:eff = pc^+i?“( l+«l^+«2 ^ 


as 


Hnr 

pc2 


a4 


jjni 

fXC^ 


(2.7) 


where the coefficient 0:4 parametrizes a possible 4PN- 
level contribution to the energy map /. 

The energy map / is determined by the requirement 
that Eq. (2.7) correctly relates the real Hamiltonian (2.2), 


2PN level, it was found that [TO] 


(|2.3|) to the EOB mass-shell condition (|2.4|), (2.5). At the 

( 2 . 8 ) 


Q!i = 2 ’ a2 = 0. 

At the 3PN level, it was further found that 
as = 0 . 


(2.9) 


One of the new results of the present 4PN-level work has 
been to find that 


0:4 = 0. 


( 2 . 10 ) 


This means that the simple energy map found at lower 
PN levels remains valid at the 4PN level. This energy 
map can be written as 


^eff 


HeS 

pLC^ 


n2 „4 


77^ — nil ~ ^2 d 
2 mi m 2 


or as 


// = Afo^A+ 2 .((( 5 -l 


III. STRATEGY 


( 2 . 11 ) 


( 2 . 12 ) 


After having recalled the basic building blocks of the 
EOB formalism, let us outline the strategy we shall use 
to transcribe the 4PN Hamiltonian of Ref. [T] within the 
EOB formalism. The need for a special strategy arises 
from the main new conceptual feature of the 4PN-level 
Hamiltonian: we found in Ref. [Tj that, contrary to pre¬ 
vious PN levels which led to ordinary (instantaneous) 
Hamiltonians, the 4PN conservative dynamics involves 
both local-in-time interaction terms (described by a usual 
Hamiltonian) and a specific (time-symmetric) nonlocal- 
in-time interaction. 

To explicitly describe this structure, it is convenient 
to henceforth replace the center-of-mass-frame variables 


Pi = —P 2 and X 12 = Xi — X 2 used in Eqs. (2.21 and (2.3) 
above by the following rescaled variables: 


X 12 Pi _ P2 

GM’ p ~ p' 


(3.1) 


Note that, in terms of these variables, the Newtonian- 
level, reduced Hamiltonian takes the simplified form 


77g-(r,p) = ip2-i. 

z r 


(3.2) 
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The (nonrelativistic, reduced) 4PN-accurate Hamilto¬ 
nian can be decomposed in two parts: 

H”[r, p] = p; s) + p; s], (3.3) 


where the first part is local in time while the second part 
is nonlocal in time. [We use brackets, iJ[r,p], instead 
of round parentheses, H{r,p), when we wish to empha¬ 
size a nonlocal functional dependence on the phase-space 
variables.] More precisely, the Hamiltonian H”Jj(r,p; s) 
is a function of phase-space variables (r, p) of the form 

^reai(r, P; s) = p) -h F(r, P) In (3.4) 


where the Hamiltonian .H< “pn is defined in Sec. V of [T] 
and the function F is defined in Eq. (3.8) there; s is a 
reduced scale with dimension of 1 /velocity^, 


._ '^phys 

■ GM ’ 


(3.5) 


where Sphys is a scale with dimension of a length. The 
Hamiltonian iJ“J/[r,p;s] is a functional of phase-space 
trajectories (r(t),p(f)). 


0 lyC 

p+OO 


/ i-oo 1 

+ (3.6) 


where Pf^ is a Hadamard partie finie with time scale 
T := 2sphys/c [see Eq. (4.2) in [1] for the definition] and 
where Fj denotes a third time derivative of the Newto¬ 
nian quadrupole moment Iij of the binary system. 




X) 


(3.7) 


In terms of the reduced variables (3.1 1 it reads 


= {GMY^x ( F F 6^^ ) . (3.8) 


Note that the nonlocal Hamiltonian p; s] slightly 

differs from what was defined in [ 1 ] as being the “non¬ 
local” part of the Hamiltonian. Indeed, there, the 4PN- 
level nonlocal piece of H was defined by taking as regu¬ 
larization scale in the partie hnie operation entering Eq. 

length 2 ri 2 /c instead of the 2 Sphys/c appearing 
As a consequence of this difference, the arbi¬ 
trary scale Sphys enters both parts, 77""'^ and of 

77"'' though it cancels out in the total Hamiltonian. Be¬ 
low, we shall separately transcribe each part, 77”'^ and 
Tjnrii^ into corresponding parts of the EOB formalism. 
We will check at the end that Sphys drops out of the final 
EOB results. 

The various techniques used in previous EOB works 
[13112] can be directly applied to transform the local 


( 13.61) tl 
in (|3.6 


part 77'"'^(r, p) of the 4PN-level Hamiltonian into its cor¬ 
responding 4PN-level EOB counterparts, A< 4 pj^, 77 < 4 pj 4 , 
and Q< 4 pn. 

By contrast, the matching of the nonlocal part 
7 jnrii[r , p] cannot be realized by employing canonical 
transformations of the usual, local kind. Here, there 
are two possibilities for incorporating the physics of 
7 jnrii[r , p] jp EOB description. The first possibility 
would be to introduce some nonlocality in the EOB build¬ 
ing blocks 5 ®®, Q. [For instance, this could be done sim¬ 
ply by transcribing the 4PN piece = C>(l/c®) into 

a corresponding nonlocal contribution — 0{l/c^) to 
Q.] However, there is a useful second possibility when 
focusing on the dynamics of gravitationally bound sys¬ 
tems. [We recall that, currently, the main application of 
the EOB formalism concerns the description of quasicir¬ 
cular, inspiralling motions (see, however, [17j for an ap¬ 
plication to scattering motions).] In that case, it happens 
that one can formally replace the nonlocality of 77"' by 
local contributions 7?^^ to the effective metric 5 ®®, 
together with an infinite string of local contributions 
to Q of the form 

00 

Q" = 5]Q2n(r')(p;)'”. (3.9) 

n—2 


In order to do so, several techniques can be used, and 
combined in various ways. In particular, one could dis¬ 
cuss the generalization of the usual, local canonical trans¬ 
formations, to nonlocal ones, with generating functions 
involving time integrals. However, it is possible to bypass 
the explicit use of nonlocal canonical transformations by 
a suitable combination of other techniques, as we shall 
now explain. 

Let us first note that 77 “^ being of order C>(l/c®) can 
be treated as a first-order effect on top of an ordinary 
(local) dynamics. Moreover, as we can neglect any cross 
effect between the various post-Newtonian contributions 
to 77"'i = + 0(l/c2) and 77"'" = C>(l/c®), it is 


enough to consider 77 “ ^ as a firs t-order addition to the 
Newtonian Hamiltonian, Eq. (3.2). Let us then consider, 
in some generality, a dynamics of the form 


77[r,p] = 77o(r,p) -Fe77i[r,p], 


(3.10) 


where 77o(r, p) = 5 p^ — y and where 77i[r, p] is either a 
higher-order Hamiltonian, involving a certain number of 
derivatives of the phase-space variables, or a nonlocal-in- 
time Hamiltonian involving integrals over time-separated 
phase-space variables [such as Eq. (3.6)]. One can for¬ 
mally consider that the second, fully nonlocal, case is 
equivalent to including an infinite number of time deriva¬ 
tives: 77i(r(t),p(f),... ,r(")(t),p(")(f),...). The first 
case (involving a finite number of time derivatives) has 
been abundantly treated in the literature, in particular 
within the context of the general relativistic two-body 
problem where higher-order Lagrangians and Hamiltoni¬ 
ans naturally occur beyond the IPN level |23]2S]' C)f 
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most direct relevance here is the work of Ref. [55] which 
considered the order reduction of the higher-order Hamil¬ 
tonian Hi{r, p, f, p) arising at the 3PN level in Arnowitt- 
Deser-Misner coordinates. Let us recall the crucial point 
made there. The replacement, within eiJi(r, p, r, p), of 
r and p by the equations-of-motion-related expressions 
r = p -I- dSn/Sp + 0{e), p = -r/r^ - SSu/^r + 0(e) 
[where Sb = f(pdr — Hdt) denotes the Hamiltonian 
action], yields the identity 


£iLi(r,p,r,p) 


= eHr‘^ir,p)+e 


dHi 

dr 


SSh 


dHi 5 Sh 
^ dp Sr 




(3.11) 


where p) = iLi(r, p, p, —r/r^) is the naive order- 

reduced version of Hi [using the 0 {e^) equations of mo¬ 
tion]. The extra terms in ( 3.11| ) (proportional to the 
variational derivatives of the action) immediately lead to 
the shifts 


5r = e 


dHi 


<5p = —e 


dHi 

1^' 


(3.12) 


of the phase-space variables needed to transform the orig¬ 
inal (higher-order) dynamics for r, p into a dynamics for 
the shifted variables r' = r -|- Jr, p' = p -|- (5p described 
by the ordinary Hamiltonian 

H^^-^ir', p') = iJo(r', p') + eHT\r', p') + 0{e^). (3.13) 


It is easily seen that a similar result holds for a 
higher-order Hamiltonian involving an arbitrary number 
of derivatives. Such a result can be further extended to a 
nonlocal Hamiltonian involving time-shifted phase-space 
variables, say r{t r), p(t -|- r), if one uses some inte¬ 
gral representation of r{t -\- r), p{t + r) in terms of r(t), 
p(t) and of the extra terms 6Sh/Sp, SSh/S r appearing 
in the equations of motion written above. (As usual, one 
can neglect “double zero” terms, i.e. work linearly in the 
variational derivatives of Sh-) 

Summarizing: Modulo some (nonlocal) shifts of the 
phase-space variables (which, in principle, can be read 
off from the order reduction procedure itsel f), one can 
reduce a nonlocal dynamics of the form (3.101 to an ordi¬ 
nary (local) dynamics of the form iL''®‘^(r, p) = Lfo(r, p) + 
eH];®‘^(r, p) -I- 0[e^), in which H[®'^(r, p) is the naive 
order-reduced version of i?i[r, pj. 


IV. DELAUNAY (ACTION-ANGLE) 
IMPLEMENTATION OF THE STRATEGY 


In order to apply this strategy, one needs, however, 
an explicit way of solving the zeroth-order equations of 
motion [so as to compute r{t + r), p(< -|- r) in terms 
of r(t), p(t)] and of then computing the orde r-red uced 
value of the nonlocal Hamiltonian Eq. (3.6). As 

in our case the zeroth-order [0(e°)] equations of motion 


are the Newtonian equations of motion, it is convenient to 
use the (Delaunay) action-angle form of the Newtonian- 
level motion. Let us recall it (using essentially the no¬ 
tation of the book M)- It is enough to consider the 
planar case. In that case the action-angle variables are 
{C,£]Q,g). Here, the action variables £, Q are related to 
the usual Keplerian variables a (semimajor axis) and e 
(eccentricity) via 

C-—y/a, Q := \/a(l - e^), (4.1) 

or 



Note that we work here and below with the rescaled 


variables (3.1). In particular, a denotes the rescaled semi¬ 


major axis a := aphys/{GM), and the time variable corre¬ 
sponding to the variables (£, £; Q, g), and the correspond¬ 
ing Hamiltonian, 


Hnr _ ^ „2 1 _ 1 

° 2^ r 2 £ 2 ’ 


(4.3) 


is i := tphys/(GM). 

C is conjugate to the “mean anomaly” i, while Q is 
conjugate to the argument of the periastron g = to. The 
explicit expressions of the Cartesian coordinates (x, y) of 
a Newtonian motion in terms of action-angle variables 
are given by 


x{C,£]G,g) = xocosg - yosing, (4.4a) 

y{£^,£]G,g)=xosing + yoCOsg, (4.4b) 


xo = a(cosu — e), (4.4c) 


j/o = a\/1 — e2 sinu, (4.4d) 

where the “eccentric anomaly” u is the function of £ and 
e dehned by solving Kepler’s equation, 

u — esinu = £. (4.5) 


In the equations above, a and e are supposed to be ex¬ 


pressed in terms of £ and G using Eq. (4.2). The solution 


of Kepler’s equation Eq. (4.5), i.e. the explicit expression 


of u in terms of £, £ and G ■, can be written in terms of 
Bessel functions: 


2 


Jn{ne) sin(n£). 


(4.6) 


Note also the following Bessel-Fourier expansions of cos u 
and sinu [which directly enter (xq, yo) and thereby (x, y)] 


LXJ ^ 

cosu = -- -f -[J„_i(ne) - J„+i(ne)] cos(n£), 
Z n 

n—l 


(4.7a) 
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OO ^ 

sinu= > — [J„_i(ne) + J„+i(ne)] sin(n^). (4.7b) 

-n 


For completeness, we also recall the expressions involving 
the “true anomaly” / (polar angle from the periastron) 
and the radius vector r: 


Xo 


yo 


ecosTt) = 

a(l — e^) 

(4.8a) 

1 -1- e cos / ’ 

= cos / = 

COS u — e 

(4.8b) 

1 — e cos u ’ 

= sin / = 

Vl — sinw 

(4.8c) 

1 — e cos u 


The above expressions allow one to easily evaluate the 
expansions of x{i), and theref rom the components 
of the quadrupole tensor Iij{i), Eq. (13.81), as power series 
in e, and Fourier series in i. One could also have used 
the known Bessel-Fourier expansions of the components 
/y (^) [28H30]. 

Let us then consider the expression 


T{t,T) := hj{t)iijit + T), 


(4.9) 


which enters the nonlocal-in-time piece (3.61 of the 
Hamiltonian. In order to evaluate the order-reduced 
value of T(t,T) we need to use the equations of motion, 
both for computing the third time derivatives of , and 
for expressing the phase-space variables at time t -I- t in 
terms of the phase-space variables at time t. This is quite 
easy to do in action-angle variables because the zeroth- 
order equations of motion following from the Hamiltonian 
Ho{C) = —1/(2£^) are simply 


= 0, (4.10a) 


d£ 

dHo 

1 

n{c), 


dHo 

di 

dC 


dt 

dG 

dC 

dHo 

= 0, 

d^ _ 

dHo 

= 0. 

dt 

d£ 

dt 

dg 


(4.10b) 


Here, we recall that t = t/{GM) is a rescaled time. We 
have introduced the notation H(£) = for the corre¬ 
spondingly rescaled (t-time) Newtonian (anomalistic) or¬ 
bital frequency: H = GMHphys (it satisfies the rescaled 
Kepler law: ft = The fact that g, L and Q are 

constant, and that £ varies linearly with time makes it 
easy to compute Iij{t + r) in terms of the values of 
{£,g,C,Q) at time t. Namely it suffices to use (denot¬ 
ing by a prime the values at time t' = t + t) 


i' = i{t + t) = £{t) + n{C)i 


(4.11) 


where r = t/{GM), together with g' = g, C = £, and 
G' = Q. Finally, the order-reduced value of is 

given by (using d/dt = Hd/d£) 


d^U 


^ ' ’ \GMJ d£3 d£3 




(4.12) 


Inserting the expansion of Iij{£) in powers of e and in 
trigonometric functions of £ and g yields in the form 
of a series of monomials of the type 


J-(£,f) = 


/ > ~nin2n^ 
ni,n2,±n3 


e"Wos(n 2 £±n 3 Hr), (4.13) 


where ni, 712 , ns are natural integers. (Because of rota¬ 
tional invariance, and of the result g' = g, there is no 
dependence of £F on g.) 


All the terms in the expansion (4.13) containing a 


nonzero value of n 2 will, after integrating over r with 
the measure dT/|f| as indicated in Eq. (3.6), generate 


a corresponding contribution to which varies with 

£ proportionally to cos(n 2 £). At this stage, we appeal 
to the usual Delaunay technique: any term of the type 
A{C) cos(n£) in a first-order Hamiltonian perturbation 
eHi(C,£) can be eliminated by a canonical transforma¬ 
tion with generating function of the type eg{L,£) = 
eB{C)sin{n£). Indeed, 


5gH^ = {Ho{C),g} = - 


OHojC) dg 
dC d£ 


= —n H(£) B{C) cos(n£), (4.14) 

so that the choice B = A/{nGt) eliminates the term 
Acos(n£) in Hi. 

This shows that all the periodically varying terms 


(with 712 7 ^ 0) in Eq. (4.13), can be eliminated by 
a canonical transformation. This proves that one can 
finally further simplify the (order-reduced) second (non¬ 
local) part of the 4PN Hamiltonian by replacing it by its 
f-averaged value. 




p27r 


d£77reai (4.15) 


i.e., 


- 1 Hr 

(4.16) 

where F denotes the ^-average of F{£^t) [which is sim¬ 
ply obtained by dropping all the terms with 712 7 ^ 0 in 
the expansion (4.13)]. This procedure yields an averaged 
Hamiltonian which depends only on C, G (and s), 

and which is given as an expansion in powers of e. Be¬ 
cause of the averaging the latter expansion contains only 
even powers of e. 

The final step of our strategy will be to match the latter 
Hamiltonian to a corresponding piece in the EOB Hamil¬ 
tonian. This matching is naturally done by performing an 
analog Delaunay reduction of the corresponding piece in 
the EOB Hamiltonian, say To start with, H^q^ 

is given, from the mass-shell condition ( |2.4| ), ( |2.5| ), as a 
function of p '^, 1 jr' and , and contains an a priori in¬ 
finite string of powers ofp'^: Qi{r') p'^ + QQ{r') p'^ + ■ ■ ■ ■ 
As p'^ is of order of the eccentricity, and as we noticed 
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that contains only even powers of e, we see that 

(as anticipated) it is enough to include in Q only terms 


even in p^. 

Having explained beforehand our strategy, we shall 
successively implement its various steps in the following 
sections. 


V. SPLIT OF THE EFFECTIVE EOB HAMILTONIAN 


Henceforth we shall set c = 1 for simplicity. Let us recall that the reduced effective-one-body 4PN-accurate 
Hamiltonian (expressed in the coordinates r', p' of the effective dynamics), i.e. the solution of the EOB mass-shell 
condition, takes the following explicit form [in which D := Q := Q/p?]: 


Hes{v',p') := ’ P ^ = jA{r') ( 1 + p'^ -f (^A{r')D{r') - 1) (n' • p')^ + Q(r',p') ). 




(5.1) 


To the (local versus nonlocal) split Eq. ( |3.3| ) of the two- 
body Hamiltonian, there corresponds a (4PN-accurate) 
split of the various building blocks A, D, and Q entering 
the effective Hamiltonian of the form 


A{r') = A\r')+A^\r') 


^>0 = 1-3 + 33+1^ 


04 4.c + «iinlnr' 


r/5 


40 +>in + 


A^\r') = 

r 

D{r') = D\r') A D^\r'), 

^V) = 1 + !S + ^| + ' 


4 in Inr' 


D^\r') = 


+ d}\ In r' 


r'4 


-‘4,111 ' 


(5.2a) 
, (5.2b) 

(5.2c) 

(5.2d) 

(5.2e) 

(5.2f) 


In Eqs. (5.2) we have used the explicit values of the 2PN- 
accurate coefficients in the potentials A} and the two 
coefficients 04 and dj, are at 3PN level, and the remaining 
ejght coefficients (o^c o^in’ « 0 in’ 41,c> 4ij„, dl\, 
level. 

The corresponding split of the “non geod esic” term Q 
of the effective-one-body Hamiltonian (5.1) is taken with 


the following structure (which will be checked to be ad¬ 
equate) : 


Q[t', p') = gOr', P') + P'). 

442 , lia.c+ 943,111 In r' 


Q (rOpO = 


r ,?2 


,,/3 


(5.3a) 

(n+pT 


962,c+ 962,111 In + 

-(n •p)^ 


r'- 


(5.3b) 


SII/ 1 11 _ 943,c + 943,111 In r' , 


O^ 0 r',p') = 


r-lZ 


w-p'y 


+ 


94,c + 960inlnr' 

-32-(n - P) 


+ O((n+p04. 


(5.3c) 


Here the coefficient (742 represents the 3PN order, while 
the eight coefficients 94,110 962,0 94,iio 

g4,o nnd in nre the 4PN level. 

In the formulas above, all the 3PN-level coefficients 
have been determined in our previous work | 12 j , namely, 


= 


94 4l7r2 


3 32 

J 3 = b2v — 
g42 = 2(4 — iv)v. 


When inserting the above I -f H split of all the func¬ 
tions A, D, and Q entering the effective Hamiltonian, one 
obtains, after expanding the right-hand side of Eq. (5.1) 


into a Taylor series with respect to 1/c^ p'^ ~ 1/r'), 

the following split of HeS- 

iLeff(r', p') = >ff(r', p') + H4(r', p') + 0(5PN). (5.4) 

Here the first effective Hamiltonian Hlg is computed from 
the functions D^, and (as if A^^, and were 
equal to zero), while the second effective Hamiltonian 
Hlff is linear in the 4PN-level functions A^^, and 
Note that both parts of the effective Hamiltonian 
are linear in the 4PN-level contributions to the functions 
A, D, and Q. For instance, the (simpler) second part of 
the effective Hamiltonian reads (at the 4PN accuracy) 


Hlh 


(E, pO = i [A^\r') + D^\r'){W ■ p')^ + g"(r', p')) . 

(5.5) 
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VI. MATCHING OF THE FIRST (LOCAL) 
PART OF THE HAMILTONIAN TO THE FIRST 
PART OF THE EOB DYNAMICS 


In this section, we shall consider the matching of the 
first part of the two- body Hamiltonian, i.e. the (local) 
part in Eq. (3.3), to the corresponding first part of 
the EOB dynamics, described by the functions Z)\ 
and discussed in the previous section. (Note that 
this first part involves all the known, lower PN contri¬ 
butions.) This matching could be done by any of the 
various techniques used in previous EOB works [Bill]. 
Here, we shall use the technique of [T2] based on requiring 
that the EOB phase-space variables (r', p') differ from 
the original (rescaled) ones (r, p) by a canonical trans¬ 
formation. It is at this stage that one should a priori 
allow for a general energy map Eq. (2.7), possibly in¬ 
volving a new, 4PN-level parameter a^, between the real 
two-body energy and the effective one. In other words, 
the two-body/EOB matching should a priori be done by 
writing that 


(i7eff(r', p'))^ ^ P) (l + P) 

+ a3(-ffr"al(l’, P))^ + a4(i7"al(l’, P))^)) > (^-l) 

and by splitting such a matching equation in two parts, 
according to our general I -I- H split. However, we already 


mentioned above that one of our results at the 4PN level 
is that the energy map Eq. (2.7) between the real two- 
body energy and the effective one does not need to be 
changed compared to previous E OB re sults [inilll], i.e. 
that we have simply 04 = 0, Eq. (2.10). 


Matching the part of the two-body Hamiltonian 

to the part H\q of the effective EOB Hamiltonian means 
that the equality 


p'))^ = p) (1 + P))^ 


( 6 . 2 ) 


holds modulo a canonical transformation between the 
phase-space variables (r,p), (r',p') given by some un¬ 
known generating function with the symbolic structure 


P') = G<3PN(r, p') + (r • p')(l + In?’) 

X (^(p'2)4 + i ((p'2)3 + ^ (6.3) 

where Gopn is known |12j . and where the omitted coef¬ 
ficients entering the 4PN level are to be found. 

We found a unique solution to this 4PN-level matching. 
Let us only give here the results for the physically most 
relevant information, i.e. the first parts of the various 
EOB building blocks. They are 


^5.c 


4237 . 22757r2 128 

T 


60 


512 


di. = 


/7243 2376l7r2 1184. 


V 45 


1536 


15 


In s I V + 


4l7r^ 

~w ~ 

/1237r2 

I 16 


221 


T 128 

as.in = 


- 260 di,i„ = 


1184 


9 L.C = 20?/ - 83?/^ -f 10?/^, g43.i„ = 0, (?62.c =-\v, q62.in = 0. 


(6.4a) 

(6.4b) 

(6.4c) 


VH. MATCHING OF THE SECOND 
(NONLOCAL) PART OF THE HAMILTONIAN 
TO THE SECOND PART OF THE EOB 
DYNAMICS 


We have explained above our strategy for the more sub¬ 
tle nonlocal-in-time second part of the two-body Hamil¬ 
tonian. It involves three steps: (i) one must explicitly 
compute the Gaveraged, o rder- reduced nonlocal part of 
the two-body Hamiltonian (3.6); (ii) one must separately 
compute the Gaverage of the second part of the ef¬ 
fective EOB Hamiltonian; and then (iii) one must match 
the two results. 


We have indicated above the Bessel-expansion tech¬ 


nique allowing one to compute to any preassigned 

order in eccentricity. This technique can be applied in 
vario us ways, e.g.: (i) by inserting the e expansion of Eq. 
(4.6), i.e. (up to e® for illustration) 


= £- 


1 3 
® 8 ® 


1 

1 ^ 


sin^ 


- 1 sin 

2 6 


3 3 27 

- - e' 

8 128 


125 

+ sin5£ + G(e®), 


1 


® sin3€+ -e'‘sin4£ 


(7.1) 


















directly into Eqs. (4.4a)-(4.4d) to compute the e- 


expansio n of t he quad rupole moment (€, e); or (ii) to 
use Eqs. (4.7a)- p?7b ) as intermediate results; or (iii) to 
use the Bessel-Fourier expansion of the quadrupole mo¬ 
ment [25^50] . say (where e = 2.718281828... should be 
distinguished from the eccentricity e) 


+ 00 




(7.2) 


We used, as a check on the results, several of these tech¬ 
niques, and we have pu shed the calculation to the seventh 
order in e. Using Eq. (7.21, we see th at the structure of 
the Uaverage ^ of J^{i, f), Eq. (4.12), will be 


grals (here written in their one-sided forms) of the type 
[see, e.g., Eq. (5.8) in [1]], 

r+oo 

Pir — cos(wr) =-( 7 e- f ln(|w|T)), (7.4) 

Jo 

where w = pPl is some integer multiple of U. This yields 
the series representation 


Keltic,g;s) 


4 / U y 




(7.5) 


Hr) = 


f PI 


\GM 


i +00 

P— — 00 


(Similarly looking, but different, series appear in tail con¬ 
tributions to gravitational wave fluxes [30].) Let us define 
the dimensionless reduced quadrupole moment 


= 2 


/ n 


6 oo 


(^) H^'^l^b(p)Pcos(pUT). (7.3) 

^ ' p—1 


The computation of the averaged Hamiltonian 
given by Eq. (4.16), is then reduced to a series of inte- 


Then, making use of the Bessel-Fourier expansion of the 
individual components of the quadrupole moment [see, 
e.g., Eqs. (A3) in [3D]], one can show that 


|4(7')P = (1 “ (1 +4 -4e^) 

- 3e(l - e^)(2(p -h 1)^ - p{2p + 3)e^) Jp_i(ep) Jp(ep) 

-b (6(p -b if - 3{p + l){5p + 2)e^ + (I2p2 + + i) - 3p^e^) Jp{ep)'^. (7.7) 

Because of the factor in the denominator, this (simplified) expression hides the fact that the term of index p is of 
order G(e^lP“^l). To see this, one can use the following equivalent (but more complicated) expression: 

|4(P)P = “ e^)(5e(jp+i(ep) - Jp_i(ep)) -b 2(6^ -b l)(jp-2(ep) - Jp+2(ep)) -be(jp+3(ep) - Jp-fep))^ 

-b(4(9e^-bl)Jp(ep)-3e(3e^-b7)(jp_i(ep)-bJp+i(ep))-b6(e^-bl)(Jp-2(ep)-bJp+2(ep))-be(e^-3)(jp_3(ep)-bJp+3(ep))^ 

-b ^4(6e^ - 1) Jp(ep) - 3e(2e^ -b 3)( Jp_i(ep) -b Jp+i(ep)) -b 6( Jp_2(ep) -b Jp+2(ep)) -b e(2e^ - 3) {jp-fep) + Jp+fep))^ 

-b (-4(3e^-b2)Jp(ep)-b3e(e^-b4)(jp_i(ep)-b Jp+i(ep)) - 6e^(jp_2(ep)-b Jp+2(ep)) -b e^( Jp_3(ep)-b Jp+3(ep))) |. 

(7.8) 


The final result for the averaged Hamiltonian reads 


i7"ai (/:,^;s) = 


real 

^ ^64 „ , / e'^^s 

y f21n2 + ln f ^ 




296 

91 — In2-b 7291n3- 


/29966 


In 2 — 


13851 

20 


In 3-b 242 In 


/ e 


,7b , 


\ cC^ 


1256 , , 2 


116722, „ 419661, „ 1953125, . 1526, /'e'>^s 


15 


In 2 


320 


ln3 - 


576 


ln5 -b 


3 


e®+G(e«) . 


(7.9) 
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We recall that the eccentricity e is considered to be the 
function of £ and Q given by Eq. (4.2). 

Let us now consider the second step, i.e. the ^-average 
of the part Hlg of the Taylor-expanded Hamiltonian HeS, 
expressed in terms of Delaunay variables C, G, i and 
g. In view of the explicit expression of Hlg (with the 
above-given forms of and Q^^), the result de¬ 


pends on the ^average of monomials involving various 
powers of 1/r' and pj, (and sometimes a logarithm of r'). 
All those computations (which must be done along the 
Newtonian motion) can be straightforwardly performed 
by expanding all formulas in the eccentricity e up to the 
order e^, using the standard Keplerian-motion results re¬ 
called above. We found 


^271- 




2/:io 


+ a"in ln(£^) -f - (20^^^ - -f 2dl]^ + 2 ( 104 ^ 1 ^ -f ln(/:2)) 


+ + 15^4 c + 3(?43,c^ + g -f g43_in) ln(£^)^ 

+ (ife (5040a”c - 5018aiii„ + 840dl^, - 533dll, + + 60qH^,) 

+ ^ (420a“ „ + 70dl\^ + + 0(e«)^ . (7.10) 

The last step is to match the contributions of these two averaged dynamics to the real (rather than effective) 
Hamiltonian. Here, we do not need to invoke any canonical transformation because the action variables of the real 
and effective dynamics are to be identified nni- One must still take into account the a priori nontrivial energy map 
relating the real and effective energies. However, it is easily seen that, as we are here dealing, on both sides, with an 
additional contribution of the 4PN level, it is enough to impose the requirement 




The unique result of this matching then gives 


TT 128, ^ TT 192 

“s.c =-^( 7 e + 21n2-f lns)i/, agjn =- —v, 

-TT / 864 1184 6496, „ 2916, 1184, 

I 5 15 ' 15 5 15 


-II 592 
7 a4.1n =- 


94^0 = 


<l62,c = 


5608 496256, „ 33048 

-^-ln2- 

15 45 5 


ln3^z/, g“ 3 i„ = 0. 


4108 2358912, 1399437, 390625 

-ln2H-ln3H- 

15 25 50 18 


In 5^^, g6kin = 0- 


(7.11) 

(7.12a) 

(7.12b) 

(7.12c) 

(7.12d) 


VIII. FINAL RESULTS FOR THE 4PN-ACCURATE EOB DYNAMICS 


Adding the results (6.4) and (7.12) of the matching procedures of the two parts of the dynamics, finally yields the 


4PN-accurate form of the EOB functions A, D, and Q. They read (we use here u= 1/r') 

'94 4l7r^' 


A{u) = 1 — 2u + 2i/u^ + 


3 32 

22757r2 4237 128 


512 


256 , „ , 

60 + 


4l7r'^ 


221 


2 64 , 

/ -I- 1 / Ini 


(8.1a) 
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D{u) = 1 + + (52j^ — 6 z/^) 

533 2376l7r2 1184 6496, „ 2916, 


/1237r2 \ „ 592 , 

( —- 260 I z/^ + -r^vlnu 


r',p')=( 


(5(r', p') = ( 2(4 - 3z/)z/ u +-— + 


5308 496256 


+ - 


827 2358912 


25 


15 


In 2 


45 


In 2 — 


33048 


ln3 — 83z^ + lOz^'^ u 


1399437 

50 


In 3 


390625 


18 


27 

ln5 1 -+ 6z/ 


( 8 . 1 b) 

)(n'p')" 

(n'-pO® + OMn'-p')®]. 


( 8 . 1 c) 


Various comments are in order. Let us first empha¬ 
size that, while the arbitrary scale s entered the separate 
pieces I and II, it has (as expected) dropped out of the 
final results. Second, we note that the 4PN-accurate re¬ 
sult for A{u) is not new, but confirms the result first 
obtained in [^. The new results with this paper concern 
the 4PN contributions to the EOB potentials D{u) and 
Q{u,p'^), except for the 4PN logarithmic contribution to 
D{u) which was previously derived in Hill mini]. For 
completeness, as the derivation used in was never pub¬ 
lished, we indicate its main steps in the following section. 
Concerning the nonlogarithmic 4PN term in D{u) that is 
linear in z/, i.e. the term of order zzit^, it was emphasized 
in that it was computable via gravitational self-force 
calculations, and a numerical estimate of its value [com¬ 
bined with some A(iz)-related contributions] had been 
previously obtained in the work of Ref. [H] (with further 
processing in Ref. [H]). For comparison, let us start by 
noticing that the numerical values of the coefficients of 
the 4PN-accurate function D{u) analytically obtained by 
us read 

D{u) = 1-1- Qvv?' + (52z/ — 6 z^^)u^ 

-k (221.572 z/ - 184.127Z/2 39.4667 z/In u)u^. 

( 8 . 2 ) 

Let us compare the analytical value of the coefficient 
of the term of order i.e. « 221.572, to the 

numerical estimate deduced from the gravitational self¬ 
force calculations of m- From the further work of [22] 
(which also made use of numerical results of [31] 1 it was 
J| num ^ 2261^, which agrees with = 221.5719912, 

within, essentially, the one sigma level. Let us also 
compare to the quantity which was directly numerically 
evaluated in [ 21 ] . namely the following combination of 
3PN- and 4PN-level coefficients of the functions A{u) and 
D{u)-. 

27 - - 9 

p\ = -— -7ai + lOa^ „ - 6 ^ 3 ^,, -b dl^ + -a\ ^^ . (8.3) 

Here dl^, 05 ,,, and denote linear-in-zz subco¬ 
efficients of the coefficients d^^c, 0 , 5 ,a and aajni re¬ 
spectively. (Beware that 05 ,in denotes here the coefficient 


of w^lnw, while, above, Ogjn denoted the coefficients of 
r'“®lnr' = —u^lnu.) Our analytical results yield the 
following value for the quantity p\: 

p\ = 64.6406. (8.4) 

while the numerical value obtained in Eq. (49) of [H] was 

^cnum^69ll. (8-5) 

Again the agreement is essentially at the one sigma level, 
which is satisfactory. 

IX. CONTRIBUTIONS TO b(u) FROM HIGHER 
PN ORDERS: 5PN AND 5.5PN 

In this section, we shall complete the above 4PN- 
accurate results on the second EOB potential Divi) by 
deriving, using various techniques, both the 5PN and the 
5.5PN contributions to D{u). The logarithmic 5PN, as 
well as the 5.5PN, contributions will be derived analyti¬ 
cally (and exactly, as functions of v) using effective action 
techniques. By contrast, only the linear-in-zz piece in the 
nonlogarithmic 5PN contribution will be estimated below 
from the numerical gravitational self-force computations 
of precession effects of Ref. [2T] . We shall then compare 
this PN knowledge of the linear-in-zz piece in D{u) to its 
numerical determination in the strong-field domain [32] . 

A. 4PN and 5PN logarithmic contributions to the 
EOB dynamics 

Let us sketch in this subsection the derivation of the 
4PN and 5PN logarithmic terms in the two-body Hamil¬ 
tonian given by one of us [S] and previously used to com¬ 
pare the EOB formalism to numerical gravitational self¬ 
force computations of precession effects [5T] . This deriva¬ 
tion used a (fractionally IPN accurate) effective action 
technique which might be of interest for other applica¬ 
tions. 

Let us consider the coupling of a localized self- 
gravitating system to an external (relativistic) tidal field. 
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We wish to reach a fractional IPN accuracy in the cou¬ 
pling between the system and the external tidal field. To 
do so it is very convenient to use the Damour-Soffel-Xu 
(DSX) formalism [55H5B] . The DSX formalism allows one 
to treat the physically nonlinear effects present in IPN 
self-gravity by means of simple, linear equations. To do 
this one must use the following exponential parametriza- 

tion of the metric [331371: 


goo = - exp I 
Awi 

goi — -w, 

Qij — ^ij exp 




(9.1a) 

(9.1b) 

(9.1c) 


In this decomposition = {w,Wi) satisfy linear equa¬ 
tions, so that we can linearly decompose as = 
w'^ + w^, where describes the self-gravitational field 
(at IPN accuracy) and describes a general external 
tidal field (also treated at IPN fractional accuracy). It 
is enough to treat the coupling to the external tidal field 
at first order in the tidal field Therefore this 

coupling is described by the action 


= 




(9.2) 


Thanks to the exponential parametrization (9.1) one sees 
that ^ is (to IPN accuracy) linear in = 

(e.g. I y/g 6goo = dw/c^), so that we get 


= 




where we denoted (as usual 133137]) 


2~'00 j^i' 


'-'Oi 


(9.3) 


:= —. (9.4) 


We can then insert in the coupling (9.3), which is linear 


in the (linear) IPN-accurate, general tidal expansion 
given by Eqs. (4.15) in (Mj- The latter expansion is a lin¬ 
ear form in the external tidal moments * (t) = 

and and their time derivatives (we 

use the multi-index notation of 
parts, one finds that the coupling (9.3) yields (modulo a 


,I37])- Operating by 


boundary term arising from a total time derivative) the 
remarkably simple result 


-S'int - X! £! / 




[yit)] 




Hr{t)Sl^[y{t)] , (9.5) 


where M™[y(t)], Sf^[y{t)] are the IPN-accurate 

Blanchet-Damour (mass and spin) multipole moments 
of the considered localized system [symbolized as y{t)] 


whose coupling to an external (tidally expanded) field we 
are considering. We recall that and S'™ are certain 
(compact-support) spatial integrals in volvi ng various mo¬ 
ments of the densities cr and cr*, Eqs. (9.4) (see 133 [37]). 
(At Newtonian order the quadrupole A/™ reduces to 
the usual quadrupole m ome nt, denoted lab above.) Note 
that the general action (9.5) can be useful for various ap¬ 


plications. For instance, it can be used to describe IPN- 
accurate radiation-reaction effects. Indeed, Refs. [3313S] 
(generalizing Ref. [40]) have shown that radiation reac¬ 
tion can be described by external-like potentials tc™, 
which are parametrized by tidal moments given by time 
derivatives of the IPN-accurate Blanchet-Damour mo¬ 
ments S™, when working at the IPN fractional 

accuracy [33] j namely 




G^eit) = 
H^it) = 


2 G bd(5) 

5 ^5 at- ’ 

A cBD(5) 

15 c5 ’ 


(9.6a) 

(9.6b) 

(9.6c) 


where the superscript, say in denotes the nth 

time derivative. One can check that the coupling defined 
by (9.5), (9.6) implies [when varying only the system vari¬ 
ables y{t) in (9.5)] the IPN-accurate radiation-reaction 

equations of motion of Ref. m_ _ 

Moreover, the same result (9.5), (]9.6[) can be used 


to compute the conservative logarithmic contributions 
arising at 4PN and 5PN. Indeed, Ref. [5] has shown 
that the first logarithmic contribution arose at 4PN 
from the fact that the (near zone) radiation reaction 
quadrupole G™(t) was the fifth derivative of a tail- 
modified quadrupole moment involving an integral 

over the past of the type (for € = 2) 


M™(t) = Mz.(t) + ^ drMf A - t) 


In 


cr 

2rn 




(9.7) 


[and similarly for 5'™(t)], where M denotes the total, 
Arnowitt-Deser-Misner mass of the local system (hence¬ 
forth treated as a nondynamical constant), and where 
ro is the length scale entering the Blanchet-Damour 
multipolar-post-Minkowskian formalism [33 ■ Note that 
the coefficient of the integral (tail) term is 4GA4/c^, i.e. 
twice the coefficient of the tail term entering the (wave 
zone) radiative multipole moments [371[33[11- The rea¬ 
son for this factor 2 is clearly general [33 so that the 
structure (9.7) applies for all values of £ (and for both 
types of moments). As a consequence, the ln(l/c) con¬ 
tribution to is given by 


M™(t) = ML{t) - 4 




Mi'A)- 


(9.8) 
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The corresponding In - contribution to the equations 
of motion is found to be conservative (though it comes 
from a tail modification of a dissipative effect) [5]. As a 
consequence of the fej/inear structure of the action (9.5), 


the 4PN and 5PN conservative logarithmic effects can 
be very simply obtained by: (i) inserting the logarith¬ 
mic contribution of (9.8) in (9.6) and then in (9.5) (with 


grext ^ ^RRin^ ^ext ^ Operating by parts to 

equalize the number of derivatives on Ml and Sl', and, 
last but not least, (iii) by adding an overall factor ^ [be- 
ca use o nly Mi[y(t)], SL[y{t)] were supposed to be varied 
in ( |9.5[ ), while we now wish to have an action where one 
has to vary two occurrences of Ml in a bilinear form 
'^Gl[Ml]Ml]- This yields the simple action 


= =-2^ Inf- 


dtJ-i;™[y(t)], (9.9) 


where 


Note that the ^-averages (here indicated by (.. .)e) have 
to be done at the IPN accuracy (which is most easily 
done by using the IPN quasi-Keplerian parametrizations 
of [IS]). After straightforward calculations, starting from 
the expression of given in |53|, one finds that the 

4PN and 5PN logarithmic contributions to A{u), D{u) 
and Q{u,pr) are (using the replacement In ^ ^ In ^ 

In u; as in 0) 


S' 

? 

II 

64 ( 

7004 

— vu Inu -|- 

- V 


5 V 

105 

D'"(m) = 

592 4 , 

/ 1420 

- vu Inw-I- 

-— 1 


15 

V 7 


/5428 592 


V 105 ^ 5 

v^ \ Mp; 


v - v 


15 


vP Inrt, 
(9.14a) 

^ Inu, 
(9.14b) 

(9.14c) 




IPN 




BD(4)x2 


he 


16 

4^ 




(9.10) 


denotes the fractionally IPN-accurate instantaneous en¬ 
ergy flux radiated as gravitational waves by the sys¬ 
tem [considered as a function of the dynamical variables 
y ^ (fa, Pa) describing the localized system, say a two- 
body system]. The corresponding Hamiltonian is simply 

Trine _ I 2GA4 /1\ ('oii'i 

'W 4 PN- 1 - 5 PN — ^ I c / 


The 4PN contributions (~ In u in A and ~ In it in 
D) agree with our complete 4PN results above, while the 
other terms are additional 5PN contributions to the EOB 
formalism. At the time of their derivation (see 0) a com¬ 
parison with the fractionally IPN-accurate computation 
by Blanchet et al. |3] of the energetics of circular binaries 
[entirely encoded in the A{u) function] showed agreement 
with A^"(ii) above for the terms linear in p, but disagree¬ 
ment for the coefficient of v^u^Xnu. However, the later 
correction, mentioned in m, to the results of | 1 | yielded 
perfect agreement with A^''(i() above [as independently 
derived in |52|, which also deduced part of the D'"(it) 
above from the combined results published in mi 


This remarkably simple result generalizes the 4PN loga¬ 
rithmic contribution of the Hamiltonian of Ref. [T] [seen 
from the near zone, the Inr term is a ~ ln(Hr/c) ~ 1 r (j) 
contribution]. It suggests that the nonlocal part 
Eq. ( |3.6| ), can be simply generalized to 

- ^ ]k .rf “ (M + U. (9.12) 


where t') is obtained from (9.10) by split¬ 
ting each squared multipole according to —>• 


(«+i) 




M 

The result ( |9.11| ) yields also a straightforward way 
to compute the logarithmic contributions to the EOB 
Hamiltonian at the 4PN and 5PN levels. A simple way 
to do it is to use the strategy we used in the text for 
namely to compute separately the ( t?-average) Delaunay 
Hamiltonian corresponding to (9.11), and the Aaveraged 
squared effective Hamiltonian, and to relate their loga¬ 
rithmic contributions via 


^In 



4 G 


V CP 


In 



.ffeff 


H 


((mi -I- m2)c^ 


(9.13) 


B. 5.5PN contributions to the EOB dynamics 

It has been recently realized that there existed half- 
integral-order PN contributions to the conservative dy¬ 
namics of binary systems As explained in m. 

such terms (which seem to violate the usual PN lore that 
conservative effects are associated with even powers of 
1 /c) have the same conceptual origin as the logarithmic 
5PN contributions discussed above; namely the fact that 
the tail-transported hereditary contribution to the near¬ 
zone gravitational field | 2 | is time dissymmetric, with¬ 
out being time antisymmetric. We shall use the same 
effective-action approach as in the previous subsection 
but extend it (as in m) to the second-order tail con¬ 
tribution. The relevant term in the near-zone metric is 
[see Eq. (79) in [47] . or, after time symmetrization, Eq. 
(5.13a) in [i^ ] 



2G 

5c^ 


x^xMB 



2 


X 



r) In 



(9.15) 
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where 


B = - 


107 

Im’ 


and where rg = rg exp(—11/12) is a rescaled ver¬ 
sion of the length scale entering the multipolar-post- 
Minkowskian formalism of Ref. [12] • (Actually, the pre¬ 
cise value of Tg is unimportant, as the scale dependence 
will disappear after taking the time symmetrization.) 

Starting from Eq. (9.15), one must time symmetrize it. 


and then insert it into the effective action [see Eq. (9.3)] 


= 




(9.17) 


As above, the time differentiated quadrupole moment 
entering this action must be considered, at this 
stage, as an external tidal moment. However, after the 
time symmetrization (and after a convenient integration 
by parts over the r integral) one finds that Si^t involves 
the double integral 


dt 


f+“ dr 




t) 


(9.18) 

where the r integral has been extended to the full range 
— 00 , -boo (which simply doubles the value of the inte¬ 
gral). 

When integrating seven times by parts the t integral, 
one sees that the bilinear form (9.18) is (despite first ap¬ 


pearances) symmetric under the exchange between the 
two functions Mij{t) and M®/*(t'). (One works modulo 
total time derivatives, as appropriate for an action den¬ 
sity.) As in the previous subsection, one can then obtain 
an autonomous action for the corresponding (conserva¬ 
tive) dynamics, simply by taking M®/* = M^-, and in¬ 
cluding an extra factor 

Finally, one ends up with an action which (after some 
operations by parts) can be written as 


^taiw _ _ ^ ^ 


B G (GM\ 


dt 

J —oo ^ 


10 c5 c3 y 

X 'M^f{t + T)-M^f{t-T) 
It could alternatively be expressed in terms of 

p /•+00 


(9.19) 


dtl_^ d.ln^(^)<)(t) 


MW(t + r) + MW(t-r) 


(9.20) 


(whose time symmetry is more evident), but the form 
(9.19) has the advantage of clearly exhibiting the in¬ 


dependence on the choice of any logarithmic time scale 
To = 2ro/c. 


Inserting in Eq. (9.19) the Fourier expansion (7.2) of 
the quadrupole moment [My = /y -|-0(l/c^)] leads to a 
(9.16) series of integrals of the form 


dr dx 

/ — sin(pHT) = sign(p) / — sinx = 7rsign(p), 

J —oo ^ J —oo ^ 

(9.21) 

where sign(p) is defined as being -bl ifp > 0, —1 ifp < 0, 
and 0 if p = 0. Finally the averaged second-order tail 
Hamiltonian is easily found to be given by 


2BGfGMy J d 
where Hphys = it/{GM) and 

OO 

p=i 


(9.22) 


(9.23) 


While the averaged 4PN Hamiltonian involved series 
of the type [see Eqs. (7.3) and (7.4)] 


^p®(alnp-b&)|/y(p)|^ 

p=i 


(9.24) 


we see that the second-order tail Hamiltonian involves 


the series (9.23) containing the seventh powers of p [or, 
more completely, in view of Eq. (9.211, the seventh pow¬ 


ers of the absolute val ues o f all the relat ive in tegers p 
entering the expansion (7.2)]. The series (9.23) entered 
the calculation by Arun et al. of the (first-order) tail con¬ 
tribution to the averaged gravitational-wave energy flux 
[30] . Using the definition (5.4) of Ref. [30], one has 


= p^{GMa)^ 32 ip{e), 


(9.25) 


where the normalized, dimensionless eccentricity func¬ 
tion ip{e) can be expanded in powers of with the result 

[301 [SO] 


<p(e) = 1 -b 


2335 


-b 


42955 

768 


e^-bO(e®). (9.26) 


Inserting Eq. ( |9.25[ ) into Eq. ( 9.22[ ) ( in which it 
is enough to use the approximations AI ~ M and 
GMHphys = ft — a“^/^) yields 


jjtair ^ 

— ~ C - 


<P(e), 


(9.27) 


Mc^ 2 (c^a)®-® 

where we recall that a = aphys/(GM), and where 

^ 128 „ 13696 

C=-—B, = +—(9.28) 


The averaged Hamiltonian (|9.27 ) is to be considered as a 


function of the action variables C and Q, via Eqs. (4.2). 
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Using the technique explained above, we can transcribe 
(9.27) within the EOB formalism by considering corre¬ 
sponding additional (tail-of-tail) 5.5PN contributions to 
the EOB potentials of the form 


A{u) = l — 2u+----\- Aq, 5 {v) u®'®, (9.29a) 


C. Linear-in-i/ piece in D{u): 5PN contribution and 
global fits 

In this subsection we consider the linear-in-i^ piece (or 
first-order self-force piece) in the expansion of the EOB 
potential D{u, v) in powers of v. 


D{u) = 1-1- 1 - D^, 5 {v) u®-®, 


(9.29b) D{u) = l + vd{u) + v'^d2{u)+0{v'^). (9.35) 


Qb.b{u,Pr) = q4i.5il^)PrU Pr + 0{p^). 

(9.29c) 

The matching between ( |9.27[ ) and the effective EOB 
Hamiltonian goes via (setting c = 1) 


For convenience, we have omitted to put a subscript 1 
on the 0{i') piece d{u). Before focusing on d{u), let us 
note that our results above have given access to the first 
two terms in the PN expansion of the 0{v‘^) piece d 2 {u), 
namely 


2 H 

V 


taiH 


M 


= (iij 


eff 


— (^6.5'*^’ + U^5.5' Pr + Q5.5(WjPr) ) • (9.30) 


The ^-averages on the right-hand side can be done at 
the Newtonian level. A convenient way of evaluating 
them is to transform them into /-averages, where / de¬ 
notes the true anomaly [see Eq. (4.8a)], using, for any 
dynamical function F 


d 2 {u) = —Qv? ^ - 260^ u'^ + 0 {u^) 

=-6- 184.1274^4-hC>(it®). (9.36) 

We note that both terms are negative, and will there¬ 
fore tend to reduce the value of D(u, v) when consider¬ 
ing comparable-mass systems. Similarly to the discussion 
in [32] of the behavior of the successive terms in the 
expansion of 

A{u, v) = 1 — 2u + va{u) + v^a 2 {u) + 0{v^), (9.37) 


{F)i = 




(9.31) 


As indicated in Eq. 


(9.29) the 5.5PN contribution to 


Q{u,Pr) formally contains an infinite string of powers 
of pr, starting at 0{p‘^). There is a one-to-one corre¬ 
spondence between {Ag 5 , 1)5 5 , ^44 5 , gga.s,...} and the 
infinite list of Taylor coefficients of the eccentricity ex¬ 
pansion of C (p{e) = C{1 + 'P 2 &^ + -\ -)■ At order 

e° we see that 


Ae. 5 {y)=Cv=^^^'Kv (9.32) 


in agreement with the results of Ref. m (see also Refs. 
[451145] ). We get the new result we were looking for by 
identifying the 0 {e^) contribution to the right-hand side 
of ( [9.30 ) to the corresponding terr n oc C ^>2 in Eq. (9.27) 
[using the O(e^) term [5D| in Eq. (9.26)]. This yields 


F 5.5 _ 619 

A^~W’ 


(9.33) 


(see also Refs. EH [sg for a discussion of other i/- 
expansions) one expects that the individual contributions 
d(u ), J 2 (w), ... will have formal singularities at the light 
ring (i.e. when u —?► |). Previous examples suggest that 


the formal singularity in d 2 (u) will be stronger, and of the 
opposite sign, than t hat in d{u). This suggests that the 
negative signs in Eq. (9.36) are related to, and indicative 
of, a singular behaviorHTfu) —>■ — 00 as u — 

Gravitational self-force studies of the precession of 
small mass-ratio orbits have allowed one to numerically 
compute the value of d{u) on a sample of orbits span¬ 
ning the interval 0.0125 < u < i [211132]. As we have 
seen above, these numerical data had led to a numerical 
estimate of the nonlogarithmic 4PN coefficient say d^ (de¬ 
noted d\ ^ above) which is consistent with our 4PN result. 
On the other hand, the determination of the higher-order 
PN coefficients, and notably of the nonlogarithmic 5PN 
coefficient dg, was extremely coarse. Namely, the tran¬ 
scription [22] of the result —6000 < < —440 [Eq. (49) 

in [ 21 ]] yielded 


- 1849 <dl< -249. (9.38) 


i.e. 

264932 

D 5 5 M = TT V = 528.4497936 v. (9.34) 

1575 

Note that both ^ 6 . 5 ( 2 /) and are proportional to 

V. This is a useful information, especially as gravitational 
self-force techniques can, currently, only give access to 
corrections linear in v. 


The additional knowledge we have acquired above 
about the analytical value of dl, as well as that of the 
5.5PN term dg.g = opens the opportunity of 

extracting more PN information from the numerical data 
on the function d{u) presented in Table VII of [32]. We 
explored various PN-like fits of the latter data on d(ii), 
keeping fixed all analytically known PN coefficients, and 
fitting for some higher-order PN expansion parameters 
involving a nonlogarithmic 5PN term d§M®, as well as a 
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combination of terms of order u® and/or u^\au and/or 
and/or u^. In order to avoid the contamination of 
higher-order PN terms (which become important in the 
strong-field domain), we have first tried to determine the 
value of from local fits, namely (error-weighted least- 
squares) fits on the weak-field interval 0.0125 < u < 0.1. 
We explored fits involving, besides dg, either one, two 
or sometimes three extra parameters. From the range of 
values of (i§ resulting from such fits, we concluded that 
the value of d§ is probably in the range 

d^num ^ _343^7g (9.39) 


This range significantly improves the previous knowledge, 
Eq. ( |9.38l ), on dg. 

We also performed some global PN-like fits of d''™(u) 

on the full available interval 0.0125 < u < h. To reach 

— — 6 

a decent goodness of fit on such a large interval, one 
needs to include at least four parameters. For instance, 
we found that a weighted least-squares fit involving 
M®, vJ, and gave a reduced equal to 2.14 for the 
following coefficients: 


dg glob = -396.905, 
dggiob = 4935.68, 
d/iob = -17634.7, 

dggiob ^ 434Q21 

In other words, the PN-like expansion 


(9.40) 


592 . 


dsi°b(y) =Qu^ + + [ 22I.57199I2-b — Inu u 


15 


-b (dg glob - 202.8571429 In m) u® 
-b 528.4497936 w®-® -b dggiob 


+ d' 


glob 7 

7 ti- 


d'gglobw®, 


(9.41) 


which combines our analytical results with the numerical 
result (9.40), yields a reasonably accurate representation 
of d{u) on the full interval 0 < m < |. In Fig. Ill we com¬ 
pare the rescaled function d (w)/(6u^) = l-b-yw-b--- 
to the corresponding numerical result d"™(M)/(6u^) 
(which was plotted in Fig. 8 of [22]). We also show, 
for comparison, some of the analytically known 3PN (i.e. 
1 -b ^ u) and 4PN corresponding results. 

Further studies are called for to see whether the new 
knowledge on the function £)(u, z/), namely [combining 
Eqs. (9.36) and ( |9.41[ )], 


D{u, z^) ~ 1 -b z^dg*°b(M) -p d 2 {u), 


(9.42) 


can be incorporated to improve the current EOB models. 
[Until now, one has only used the 3PN approximation to 
D{u,v) = 1 -b D 2 {v) -b Dziv) u®.] As was the case for 
the main EOB radial potential A(m, v), it might become 
necessary to replace the PN-like expansion of D{u^ v) by 
some resummed one. And as was the case for A{u,v)^ 



FIG. 1. Comparison of the numerical result |32| 
d“™(u)/(6«^) for the rescaled function d{u)/(Qu^) with the 
global fit d®*°®(u)/(6M^), Eq. ( |9.41[ ), and the analytically 
known 3PN |12| and 4PN results (derived above). 


it might then be necessary to appeal to specific numer¬ 
ical simulations to calibrate such a resummed form of 
D{u,v). There are several numerical experiments that 
might inform the strong-field values of D{u,i'), notably 
periastron advance studies (such as |52|) and strong-field 
scattering simulations (such as [II]). 


X. CONCLUSIONS 


Let us summarize our main results. 

We decomposed the recently derived, nonlocal-in-time 
4PN Hamilt onian [1] into two parts: (i) a local-in-time 
part Eq. ( |3.4| ), [which includes an explicit logarithmic 
contribution F'(r, p) ln(r/s), where s is a scale separating 
near-zone from wav e-zon e effects]; and (ii) a nonlocal-in¬ 
time part Eq. (3.6), (which involves an integration 


over arbitrarily large time separations). 

We transcribed the local dynamics into an equiv¬ 
alent (gauge-invariant) EOB description by performing 
a 4PN-accurate canonical transformation between the 
Arnowitt-Deser-Misner phase space and the EOB one. 
This led to EOB potentials A^(r'), &{r') and Q^(r',p(,) 
given (in exact form) by Eqs. ( |5.2b ), (5.2e), (5.3b) with 
Eqs. (6.4). Note, in particular, that is a polynomial 
in p), which involves only p'^ and p/®. 

On the other hand, we had to introduce a new tech¬ 
nique for transcribing the nonlocal part into an 
equivalent EOB description by means of ordinary (local- 
in-time) potentials A^^(r'), Z)^^(r') and (5^^(r',p(,). Our 
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technique combined several steps: (i) a formal order re¬ 
duction allowing one to express the time-separated phase- 
space variables r{t + T), p(t-|-T) as functions of r(t), p(t); 
(ii) a transformation to action-angle variables £,C]g,G', 
and, finally (and crucially) (iii) the use of Delaunay’s 
method, i.e. the elimination of periodic terms in the an¬ 
gle variables by suitable canonical transformations. Af¬ 
ter formally making an infinite number of such Delau¬ 
nay transformations, we ended up with a new Delaunay- 
like Hamiltonian depending only on the ac¬ 

tion variables C, Q (and given by an angle average of 
the original H^^). (Let us note in passing that the orig¬ 
inal derivation of the 2PN EOB dynamics m was sim¬ 
ilarly based on the 2PN-accurate Delaunay Hamiltonian 
derived in [54). 1 The resulting Delaunay Hamiltonian 
cannot be translated into explicit, polynomial 
EOB counterparts Indeed, the last poten¬ 

tial Q^^jr ',p'^) must be taken as an infinite series in p(,, 
Eq. (3.9), corresponding to the infinite expansion, Eq. 
(7.5), expressing G) in terms of the Eourier-Bessel 

expansion, Eq. (7.2), of the quadrupole moment of the 
system. Here, we have determined the value of up to 
p'^ (i.e. the maximum power of p(. entering its first part 
Q^). Our results for and are given in Eqs. 

By combining parts I and H, we could transcribe [mod¬ 
ulo Oij}'^) terms] the original nonlocal Hamiltonian into 
equivalent 4PN-accurate EOB potentials A, D and Q, 


see Eqs. (8.1). We checked that the nonlogarithmic 4PN 
contribution to A and D agree with previous (analytical 
or numerical) results. 

To complete our 4PN-level results, we used an effec¬ 
tive action technique to compute some higher-order con¬ 
tributions. More precisely, we showed how to directly 
derive the logarithmic contributions to the EOB dynam¬ 
ics arising at 4PN and 5PN, Eqs. (9.14), as well as the 
half-integral contributions first showing up at the 5.5PN 
level [Eqs. (9.32) and (9.34)]. 

Finally, we used this knowledge to improve the ex¬ 
traction of PN information from gravitational self-force 
numerical data on the linear-in-i^ piece d{u) in D{u,v). 
This led to an approximate determinatio n of t he nonlog¬ 
arithmic 5PN contribution to d(u), Eq. (9.39), and to a 
simple, PN-like representation of the behavior of d(u) in 
the strong-held domain, Eqs. dOTI and 
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